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Method and Apparatus For Reducing Three -Dimensional Shape Data 

This application is based on patent application Nos. 11-371265 
and 11-371266 filed in Japan, the contents of which are hereby 
incorporated by references . 

BACKGROUND OF THE INVENTION 

This invention relates to a method and apparatus for 
reducing three-dimensional shape data forming a polygon model. 

A huge number of measurement data are obtained by, for 
example, a three-dimensional shape measuring apparatus having 
quite a high space resolution if an image quality is high. 
Accordingly, if the data have a higher image quality than 
necessary, it disadvantageously results in a large burden on a 
storage device during succeeding processing and a considerable 
reduction in processing speed. Thus, it is necessary to ease 
this problem as much as possible by reducing the number of 
measurement data. However, in such a case, if a difference 
between a polygon model (hereinafter, referred to as 
approximated polygon model) formed by the measurement data after 
the data reduction and a polygon model formed by the actually 
measured three-dimensional shape data (hereinafter, referred to 
as original polygon model) is excessively large, the measurement 
data after the data reduction cannot be effectively utilized. 
Therefore, it is necessary to control error between the 
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approximated polygon model and the original polygon model in the 
data reduction. 

Particularly in the case, of t,hree- dimensional shape data 
obtained by a three-dimensional shape measuring apparatus, 
actually measured three-dimensional shape data are converted 
into an approximate value by the data reduction. Thus, error 
between the approximated polygon model and the original polygon 
model is preferably controlled using the actually measured 
three-dimensional shape data as a reference. 

On the other hand, various methods have been proposed to 
reduce the number of three-dimensional data forming a triangular 
polygon model by approximating a model within a desired 
permissible error range. 

For example, in "Hierarchical Approximation of Three- 
Dimensional Polygon Models" (by Horikawa and Totsuka, pp. 3-7 of 
Proceedings of the 5^*" Sony Research Forum, 1995) is disclosed a 
method for reducing the number of three-dimensional shape data 
by estimating contributions of respective edges forming a 
triangular polygon model in the form of an estimation function 
expressing a volume which changes in the case that edges are 
deleted, and deleting edges, for example, in the case that they 
are parallel to surrounding surfaces, surfaces are small and/or 
edges are short. 

Further, in "Simplification Envelopes" (by Cohen Jonathan, 
Amitabh Varshney, Dinesh Manocha, Greg Turk, Hans Weber, Pankai 
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Agarwal, Frederick Brooks and William Wright, pp. 119-128 of 
Proceedings of SIGGRAPH 96 in Computer Graphics Proceedings, 
Annual Conference Series, 1996, ACfl SIGGRAPH) is disclosed a 
method for generating models obtained by reducing an original 
polygon model by an error £ and enlarging it by the error e and 
reducing the niomber of data of the original polygon model such 
that an approximate polygon model lies between the enlarged and 
reduced models. 

Further, in "Surface Simplification Using Quadric Error 
Metrics" (by Michael Garland and Paul S. Heckbert, Carnegie 
Mellon University, Proceedings of SIGGRAPH 97 , 1997 ) is 
disclosed a method for reducing the number of three-dimensional 
shape data forming a triangular polygon model by converging two 
or more vertices (three-dimensional shape data) forming an edge 
or surface of the triangular polygon model to one vertex. In 
this method, a distance between the converged vertex due to edge 
shrinking or surface shrinking and the respective planes forming 
the original triangular polygon model influenced by the edge 
shrinking or surface shrinking is controlled as an error in the 
data reduction. The number of data is reduced within a 
specified permissible error range by performing such an edge 
shrinking or surface shrinking that this error lies within a 
specified permissible range. 

In this method , specifically , when the edge shrinking is 
applied for one edge, the error estimation is executed for the 
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other edges or surfaces of the polygon model that are to be 
influenced by the edge shrinking because the shrinking of the 
one edge changes edges and surface.? surrounding the one edge. 
After the error calculation is completed, the shrinking of the 
one edge is determined. This processing is repeated to reduce 
the number of the three-dimensional shape data within a 
predetermined permissible range. 

Since the data reducing method disclosed in "Hierarchical 
Approximation of Three -Dimensional Polygon Models" only deletes 
edges based on evaluation values thereof, an error in the 
approximation of a model by the edge deletion cannot be 
estimated although the three-dimensional data of characteristic 
portions of a polygon model remain. Thus, this method is 
effective in the case of reducing the number of data of a 
polygon model for merely displaying a three-dimensional image 
such as virtual reality and games in computer graphics . 
However, this method cannot be effective in the case that the 
error between the original measurement data and the data of the 
approximated model after the data reduction needs to be 
controlled. 

On the other hand, the data reducing method disclosed in 
"Simplification Envelopes" can strictly control such error. 
However, it is necessary to generate polygon models by enlarging 
and reducing the original polygon model by an error e and save 
the data of these polygon models during the data reduction. 
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Thus/ the number of data relating to the processing becomes 
huge, causing the problem of large burden on the memory capacity 
and calculations for the data r^ducti-on. 

Further, the method disclosed in "Surface Simplification 
Using Quadric Error Metrics" controls the distances from the 
vertices of the polygon model after the approximation from the 
surfaces forming the original polygon model as errors. Thus, 
this method is not necessarily a suitable method in the case 
that data reduction is desired to be performed by controlling 
errors of the polygon models after the approximation from the 
vertices forming the original polygon model like the measurement 
data of the three-dimensional shape measuring apparatus. 

Further, in the case of measurement data obtained by 
measuring rough surfaces or cut surfaces of a measurement object 
by means of a three-dimensional shape measuring apparatus, e.g., 
in the case that the orientation of surfaces locally suddenly 
change in a zigzag manner in a polygon model, errors for these 
zigzag surfaces are estimated. Therefore, errors are estimated 
as large values as the data reduction progresses, thereby 
presenting the problem of insufficient data reduction. 

Furthermore, each time the edge or surface shrinking is 
performed, it is necessary to estimate errors with respect to 
the edges or surfaces of the shrunk polygon model that are to be 
influenced by the edge or surface shrinking. This will increase 
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the number of calculations, consequently causing the problem of 
hindering the data reduction. 

SUMMARY OF THE INVENTION 

It is an object of the present invention to provide a 
method and apparatus for reducing three-dimensional shape data 
which are free from the problems residing in the prior art. 

It is another object of the present invention to provide a 
three-dimensional shape data reducing method and apparatus which 
enable higher speed processing. 

It is another object of the present invention is to provide 
a three-dimensional shape data reducing method and apparatus 
which enable proper data reduction based on precise error 
estimation for measurement data obtained by a three-dimensional 
shape measuring apparatus . 

According to an aspect of the invention, a method for 
reducing three-dimensional shape data, comprises calculating 
estimation values for portions to be deformed by converging two 
or more vertices of a polygon model, and reducing the nximber of 
data for the polygon model by converging two or more vertices of 
the polygon model when the estimation values are equal to or 
below the predetermined permissible value. 

The portions to be deformed may be surfaces of the polygon 
model. Estimation values for surfaces to be deformed by 
shrinking edges or surfaces may be calculated based on distances 
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between the respective surfaces and all the original vertices 
involved in the surface deformation. 

Alternatively, the number of data for the polygon model may 
be reduced by converging two or more vertices of one portion of 
the polygon model after another portion, repeatedly. Before 
each data reduction, the portion having been involved in all the 
previous data reductions is defined as a reduction prohibition 
area, and a succeeding data reduction is applied to a portion 
other than the reduction prohibition area. 

According to another aspect of the invention, an apparatus 
for reducing three-dimensional shape data, comprises an 
estimation value calculator for calculating estimation values 
for portions to be deformed by converging two or more vertices 
of a polygon model, and a data reducing device for reducing the 
number of data for the polygon model by converging two or more 
vertices of the polygon model when the estimation values are 
equal to or below the predetermined permissible value. 

The estimation value calculator may calculate estimation 
values for surfaces to be deformed by shrinking edges or 
surfaces based on distances between the respective surfaces and 
all the original vertices involved in the surface deformation. 

Also, the data reducing device may reduce the number of 
data for the polygon model by converging two or more vertices of 
one portion of the polygon model after another portion, 
repeatedly. Further, there may be provided a judging device 
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for Judging based on a calculated estimated value whether a 
portion of the polygon model is permissible for the deformation, 
and a prohibition area defining device for defining a portion 
having been involved in all the previous data reductions 
executed by the data reducing device as a prohibition area to 
keep a succeeding data reduction from being applied. 

These and other objects, features and advantages of the 
present invention will become more apparent upon a reading of 
the following detailed description and accompanying drawings. 



BRIEF DESCRIPTION OF THE DRAWINGS 

FIG. 1 is a block diagram showing a construction of a 
three-dimensional shape measurement data processing system 
according to an embodiment of the invention; 

FIG. 2 is a block diagram showing a construction of a shape 
data processor; 

FIG. 3 is a diagram showing an example of approximation of 
a polygon model by data reduction; 

FIG. 4 is a concept diagram showing a data reducing method 

by edge shrinking; 

FIG. 5 is a concept diagram showing a data reducing method 

by surface shrinking; 

FIG. 6 is a diagram showing an estimation value used in a 
data reducing method according to the embodiment; 
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FIG. 7 is a flowchart showing a procedure of a three- 
dimensional shape data reduction processing; 

FIG. 8 is a diagram showing edge shrinking in a portion 
where a discontinuous surface is formed on a polygon model by a 
random noise data; 

FIG. 9 is a diagram showing a difference between an error 
estimation value according to the embodiment and a conventional 
error estimation value in the data reduction shown in FIG. 8; 

FIG. 10 is a diagram showing shrinking of a zigzag surface 
of a polygon model having such a surface; 

FIG. 11 is a diagram showing a conventional error 
estimation value in the data reduction shown in FIG. 10; 

FIG. 12 is a diagram showing an error estimation value 
according to the embodiment in the data reduction shown in FIG. 
10; 

FIG. 13 is a flowchart showing another embodiment of the 
three-dimensional shape data reduction procedure; 

FIG. 14 is a diagram showing an example of a data reducing 
process by the operation procedure shown in FIG. 13; and 

FIG. 15 is a diagram showing a modification of the 
flowchart shown in FIG. 13. 

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS OF THE 
INVENTION 
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Referring to FIG. 1, a three-dimensional shape measurement 
data processing system (hereinafter, "shape data processing 
system") is provided with a three -.dimensional shape measuring 
apparatus 1 for obtaining three-dimensional shape data of a 
measurement object G to be measured, an external storage device 
2 for storing three-dimensional shape data obtained by the 
three-dimensional shape measuring apparatus 1, a shape data 
processor 3 for reducing the three-dimensional shape data 
obtained by the three-dimensional shape measuring apparatus 1 by 
a data reducing method to perform a specified data processing, 
an external storage device 4 for storing the three-dimensional 
shape data processed by the shape data processor 3, a display 
device 5 for displaying a three-dimensional image of the 
measurement object G such as a polygon model using the three- 
dimensional shape data after data reduction, and an external 
processing device 6 such as a CAD/CAM (computer aided design and 
computer aided manufacturing) device in which the three- 
dimensional shape data processed by the shape data processor 3 
are utilized. 

The three-dimensional shape measuring apparatus 1 measures 
a surface shape of the measurement object G using a conjugated 
autofocus system. As shown in FIG. 2, the three-dimensional 
shape measuring apparatus 1 basically includes an illuminating 
unit 101 for generating illumination light, a confocal optical 
system 102 for projecting illumination light from the 
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illuminating unit 101 onto the measurement object G and 
introducing the illumination light reflected by the measurement 
object G to a photodetector 103, and the photodetector 103 for 
receiving the illumination light reflected by the measurement 
object G and outputting it after photoelectrically converting it 
into an electrical signal, and a beam splitter 104 for splitting 
the light reflected by the measurement object G and having 
transmitted through the cpnfocal pptical system 102 
(hereinafter, referred to as merely "optical system 102"), and a 
calculator 105 for calculating coordinates of surface positions 
of the measurement object G using light reception signals 
outputted from the photodetector 103. 

The three-dimensional shape of the measurement object G is 
measured by placing the measurement object G on a measurement 
table 106 provided in front of the three-dimensional shape 
measuring apparatus 1 and measuring the surface shape of the 
measurement object G at each height position while elevating and 
lowering the measurement table 106 at specified intervals along 
height direction. Specifically, assuming that Y-, Z- and X- 
directions denote the height direction of the measurement object 
G, the optical axis direction of the optical system 102 and a 
direction normal to Y-direction and Z -direction, respectively, 
the three-dimensional shape of the measurement object G is 
measured by measuring Z -coordinate of the measurement object G 
at each Y-coordinate while changing the Y-coordinate of the 
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measurement object G at specified intervals. The three- 
dimensional shape data (x, y, z) of the measurement object G 
measured by the three-dimensional shape measuring apparatus 1 
are inputted to the shape data processor 3. 

The three-dimensional shape measuring apparatus 1 is 
provided with a photoelectrically converting device 103a 
arranged in parallel to X-axis and formed of a line sensor, and 
the coordinate ( Z -coordinate) of a surface position of the 
measurement object G . is measured using a light reception signal 
fQj^ each pixel (each pixel position corresponds to an X- 
coordinate) of the photoelectrically converting device 103a. An 
^ illumination light in the form of a slit extending in X- 

^ direction is emitted from the illuminating unit 101 so that the 

light reflected by the measurement object G is incident on the 
Ljl entire photoelectrically converting device 103a. 

'p. The Z -coordinates of the measurement object G are measured 

as follows. First, illumination light in the form of the slit 
extending in X-direction is emitted from the illuminating unit 
101 provided at an image-side focusing position of the optical 
system 102, This illumination light is gathered at an object- 
side focusing position of the optical system 102 to form a line 
beam extending in X-direction at this focusing position . 
Subsequently, the photoelectrically converting device 103a is 
repeatedly exposed in a specific cycle while the object-side 
focusing position (i.e., line beam) of the optical system 102 is 
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moved in the optical axis direction until this focusing position 
reaches the inside of the measurement object G. An exposure 
position where the level of the light reception signal is at 
maximvim is calculated by the calculator 105 for each pixel of 
the photoelectrically converting device 103a using the levels of 
the light reception signals obtained at the each exposure 
position, and this position is assumed to be the Z -coordinate of 
the outer surface of the measurement object G. 

The shape data processor 3 is provided with a polygon 
generator 31, a data reducing device 32, a permissible value 
input device 33 and a data output device 34. The polygon 
generator 31 generates a polygon model of the measurement object 
G based on the three-dimensional shape data (x, y, z) of the 
measurement object G inputted from the three-dimensional shape 
measuring apparatus 1. The data reducing device 32 reduces the 
three-dimensional shape data by edge shrinking or surface 
shrinking of the polygon model- In other words, as shown in 
FIG. 3, a precise polygon model A based on the measurement data 
is approximated to a polygon model B within a specified 
permissible error range by suitably reducing the vertices 
(three-dimensional shape data) of the polygon. 

The data reducing device 32 includes an estimation value 
calculator 321, a estimation value sorter 322 and an 
edge/surface shrinking section 323 to carry out a data reduction 
processing to be described later. 
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The estimation value calculator 321 calculates estimation 
values used to judge whether the respective edges or surfaces 
can be shrunk before shrinking edges or surfaces forming the 
polygon model by converging two or more vertices to one vertex. 
This estimation value is described later. 

The estimation value sorter 322 includes a calculator for 
calculating a minimum estimation value by rearranging a 
plurality of estimation values calculated for the respective 
edges or surfaces by the estimation value calculator 321 in an 
increasing or decreasing order. The edge/surface shrinking 
section 323 compares the minimum estimation value calculated by 
the calculator of the estimation value sorter 322 with a 
permissible value inputted via the permissible value input 
device 33 and performs an edge or surface shrinking 
corresponding to the minimum estimation value when the minimum 
estimation value is equal to or below the permissible value. 
The estimation values and an edge shrinking method are described 
later. 

The permissible value input device 33 is operated by an 
operator to input the permissible value used to judge whether 
the edges or surfaces can be shrunk as described above. This 
permissible value corresponds to a permissible range within 
which the edge or surface shrinking is assumed to have a little 
influence on characterizing portion(s) of the polygon model in 
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approximating the polygon model, and is suitably inputted by the 
operator according to the shape of the measurement object G. 

The data output device 34 is an interface for outputting 
the reduced three-dimensional shape data (three-dimensional 
shape data forming the approximated polygon data) to the 
external storage device 4. the display device 5 and the external 
processing device 6. 

At least part of devices of the shape data processor 3 may 
be provided in another independent processing unit, and their 
functions may be realized by cooperation of the processing units 
and the shape data processor 3. Alternatively, these devices 
may be processing units which are separately provided in the 
shape data processor 3 and cooperate with each other. 
Conversely, at least two functions of these devices may be 
realized by a processing program of one processor, e.g.. the 
shape data processor 3 or a processor provided therein. 

The external storage device 4 is adapted to store the data 
in an external storage medium such as a magnetic disk, an 
optical disk or a magnetooptical disk. Upon input of the three- 
dimensional shape data from the shape data processor 3, the 
external storage device 4 stores the data in an external storage 
medium in a specified file format. 

The display device 5 is an electronic display device such 
as a cathode-ray tube or a liquid crystal display device. Upon 
input of the three-dimensional shape data from the shape data 
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processor 3, an Image of the approximated polygon model of the 
measurement object G is displayed on the display device 5 based 
on the inputted data (see the p^olygon model B of FIG. 3). 

The external processing device 6 is a device for 
fabricating, for example, a mold or a model having substantially 
the same configuration as the surface configuration of the 
measurement object G. Upon input of the three-dimensional shape 
data from the shape data processor 3, a mold or the like is 
automatically fabricated by cutting a work based on the inputted 
three-dimensional shape data. 

Referring back to FIG, 1, at least part of the measuring 
apparatus 1 and the devices 2, 4, 5, 6 may be incorporated in 
the shape data processor 3. 

Referring to FIGS. 4 and 5, next, the data reducing method 
will be described. The data reducing method by edge shrinking 
is such a method for reducing the number of data that, in the 
case that a polygon model is defined by vertices Vi to Vio as 
shown in FIG. 4, two vertices Vi, V2 are, for example, so 
converged to a vertex Vi' as to shrink an edge Edi2 . connecting 
the vertices Vi, V2. 

In the data reducing method by edge shrinking (hereinafter, 
referred to as "edge shrinking method"), estimation values for 
estimating the influence of the edge shrinking on a change of 
the configuration of the polygon model for all edges Edij (i, j = 
1, 2, ...10, combinations of vertices adjacent to each other 
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when Iptj) are first calcuiated, and the vertices Vi, Vj at the 
opposite ends of the edge Edi^ corresponding to a minimum 
estimation value are converged to a specified vertex Vi' for 
data reduction if the minimum estimation value is equal to or 
below a permissible value. 

Subsequently, the edge shrinking is performed in a similar 
manner for the polygon model (approximated polygon model) in 
which the vertices Vi, Vj are converted to the vertex Vi' , 
thereby further reducing the data. Thereafter, the edge 
shrinking is performed edge by edge in a similar method until no 
more edge shrinking can be performed for the polygon model after 
edge shrinking (i.e., the approximation of the polygon model is 
repeated until the minimum estimation value becomes larger than 
the permissible value) to reduce the data. 

It should be noted that the vertex Vi' is an approximated 
point set in a specified position between the vertices Vi, v-, 
such as a middle point of the edge Edij as disclosed in the 
aforementioned paper "Hierarchical Approximation of Three- 
Dimensional Polygon Models" - 

The data reducing method by surface shrinking (hereinafter, 
referred to as merely "surface shrinking method") is such a 
method for reducing the number of data that, in the case that a 
polygon model is defined by vertices Vi to Vn as shown in FIG. 
5, vertices Vi, V2, V3 are so converged to a vertex Vi' as to 
delete a surface P123 enclosed by the vertices Vi, vj, V3. 
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Basically, the data reducing method by surface shrinking is 
carried out in a procedure similar to that of the data reducing 
method by edge shrinking. Spfecif ioally . estimation values for 
estimating the influence of the surface shrinking on a change of 
the configuration of the polygon model for all surfaces Pijk (i, 

j, k = 1. 2, 10, combinations of vertices adjacent to each 

other when i^j=*k) are first calculated, and the vertices Vi. Vj. 
Vk defining the surface P^k corresponding to a minimum 
estimation value are converged to a specified vertex Vi' for 
data reduction if the minimum estimation value is equal to or 
below a permissible value. Thereafter, the surface shrinking is 
performed surface by surface in a similar method until no more 
surface shrinking can be performed for the polygon model after 
surface shrinking to reduce the data. 

In the case that the three-dimensional shape data are 
successively reduced by the edge reducing method or surface 
reducing method, the three-dimensional shape data forming the 
measurement object G after the data reduction processing are 
different from those of the measurement object G actually 
measured by the three-dimensional shape measuring apparatus 1 
and serve as approximated three-dimensional shape data. Thus, 
it is necessary to perform the data reduction processing while 
controlling an error of the polygon model after the data 
reduction (hereinafter, referred to as "approximated polygon 
model") from the polygon model formed by the actually measured 
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three-dimensional shape data (hereinafter, referred to as 
"original polygon model" ) . 

In the data reducing method according to this embodiment , 
if it is assumed that a triangular polygon Pabc is approximated 
to a triangular polygon Pabc' by* for example, approximating 
vertices a, b, c (actually measured three-dimensional shape 
data) to vertices a', b' , c' as shown in FIG. 6 through, for 
example, repetition of the edge shrinking, distances of a 
surface including the triangular polygon Pabc' after the 
approximation from the original vertices a, b, c involved in the 
approximation of this surface are controlled as errors. 

In the example of FIG. 6, the vertices involved in the 
approximation are assumed to be the vertices a, b, c in order to 
simplify the description. However, in an actual processing, 
vertices other than the vertex a are involved in the 
approximation, for example, while the vertex a is converged to 
the vertex a'. Therefore, the above errors include distances of 
these vertices from the triangular polygon Pabc' • 

Next, the estimation value of the error used in the data 
reducing method according to this embodiment is described. 
Generally, a surface p in a three-dimensional space is expressed 
by Equation (1): 

yljc+5y + Cz + D-0 ..-(1) 
where y4^+5'^ + C' =1. 
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On the other hand, if (v^, Vy, v^) denotes the coordinates 
of point V in the three-dimensional space, a distance d between 
the point v and a surface p is calculated by Equation (2): 



d =^v, +Cv^ +Z) 



-(2) 



and can also be expressed by matrices as in Equation (3): 
d = P^ V ...(3) 



where P ~ 



'a 






b 








v = 




c 






d 




1 



±s a transposed matrix of P. 



Accordingly, the square of the distance d is expressed by 
Equation ( 4 ) : 

= (P^ -K)' = (P"^ 'V)'(V^ P) = {V'V^yP = P"^ Kv'P • . . (4) 
where is a transposed matrix of V 
Kv=V'V'^ 



V V V V 



^ x^ y ^y 

V V V V 

^ z ^ y^ z 



V V 

^ y^ z 



^x 

v., 



Equation (4) expresses the square of the distance d between 
one plane p and one point v not included in the plane p. In the 
case that a plurality of such points Vi, Vz, . . .Vm exist for the 
plane p, a sum A(p) of the squares of these points Vi (i ~ 1, 
2, ,..m) from the plane p is expressed by Equation (5): 

m m m 

A(p)=grf/' -^(^"^ iCv, •F) = P^-g/s:v-P . ..(5) 

In the data reducing method according to this embodiment, 
if there are n planes of a polygon model which are influenced by 
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the edge shrinking, for example, when data are reduced by edge 
shrinking, distances of a plurality of original vertices 
(original three-dimensional shape data) involved in the data 
reduction processing for the approximation from the planes are 
estimated for all such planes. Accordingly, the sum A{p) of the 
squares of the distances is calculated by Equation (5) for each 
of such planes, and e which is obtained by dividing a sum of the 
sums A(pi) by a total number M of all vertices involved in the 
data reduction processing as seen in Equation (6) is assumed as 
an error estimation value: 




Further, Equation (7) is obtained by substituting 

m 

F^ -^XVf Pof Equation (5) for A(/^) in Equation (6): 




I"' 

It should be noted that in Equation (7), 

P- denotes a coliamn vector for an equation of a plane pi 
subject to the deformation; 

pI denotes an inverse matrix of P^ ; 

Kv^ denotes an error matrix for a j-th vertex Vj involved 
in the deformation of the plane pi when an error matrix for data 
of an original vertex v involved in the deformation of the plane 
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Pi ^ is defined by Kv = V'V^ ( K denotes a column vector for the 
data) ; 

denotes a sum of error matrix Kvj for a vertex Vj 
involved in only the deformation of the i-th plane pi; and 
n denotes a svun of planes pi to be deformed. 

In Equation (6), A(pi) is a sum of the squares of all 
original vertices involved in the approximation of the i-th 
plane pi from the plane pi, and mi is the number of terms in the 
sum of the squares of the distances . 

In the calculation by Equation (6), since a matrix Kvi is a 
symmetric matrix as shovm in Equation (4), only diagonal 
elements and elements above them may be stored in the memory for 
the matrix Kvi calculated for the original vertex Vi. This 
enables a memory capacity necessary for the calculation to be 
reduced. Further, the respective terms mi in the denominator of 
Equation (6) are the number of terms of the sum of the squares 
of the distances in the i-th approximated plane pi. Since the 
element at the right bottom corner of the matrix Kv is 
constantly "1" and the element at the right bottom corner of a 
matrix 2Kv in the calculation of the sum of the squares of the 
distances is mi, calculation by Equation (6) may be conducted 
using the value of this element. Thus, the calculation is 
easier because it is not necessary to separately calculate the 
number of terms mi. 
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Further, it may be appreciated to use weighted error matrix 



provided a weight factor for each vertex Vi in consideration 
of the influence of the vertex Vi to an approximated polygon 
model or other causes. The weighted error matrix is obtained by 
multiplying the error matrix for the vertex Vi by a 

predetermined weight factor iv- for the vertex Vi as >v, Kv^ . 

In the case that Equation (5) is applied in estimating an 
error for a certain plane p, if original vertices Vi involved in 
the processing of edge shrinking or surface shrinking for the 
plane p are specified, a sum SKvi of the matrices Kvi obtained by 
the three-dimensional shape data of these vertices Vi serves as 
a parameter indicating a degree of the error in the edge 
shrinking or surface shrinking. In the description below, 2Kvi 
is referred to as an error matrix. 

Next, a procedure of the three-dimensional shape data 
reduction processing is described with reference to FIG. 7. The 
three-dimensional shape data of the measurement object G 
obtained by the three-dimensional shape measuring apparatus 1 
are inputted to the shape data processor 3, in which a 
triangular polygon model is generated in the polygon generator 
31 (Step #1). Further, an operator inputs a permissible value 
used to Judge whether the edge shrinking processing can be 
performed via the permissible value input device 33 (Step #3). 



for the calculation of the error e . 



Specifically, there are 
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Subsequently, estimation values (j = 1, 2, 



. . .N, N is a 



total number of edges) in the case that all edges of the 
original polygon model are shrunk are calculated by Equation (6) 
in the estimation value calculator 321 (Step #5). For instance, 
in the example of FIG. 4, estimation values e-, (j = 1, 2, ...19) 
in the case of shrinking the edges Edir (r = 2 to 5, 9, 10), Edas 
(s = 5, 6, ...9) and Edtt+i (t = 3, 4, ...10, where t+1 = 11 is 
assumed to be t+1 =3) are calculated. 

At this time, if the edge Edi2 is shrunk by, for example, 
converging the vertices Vi, Vz to the vertex Vi', the original 
planes are all deformed and eight planes pi, P2, --.Pio 
(hereinafter, planes after the deformation are referred to as 
approximated planes) are formed by the vertex Vi' and the 
vertices V3 to Vio not involved in the edge shrinking. A sum of 
the squares of the distances A(Pk) (k = 1, 2, ...8) are 
calculated for each of these approximated planes Pk, a total 
number M of the elements included in all the sums A(Pk) is 
calculated, and an evaluation value e(Edi2) in the shrinking of 
the edge Edi2 is calculated by Equation (6). The evaluation 
value e(Edi2) is expressed by Equation (8): 



At this time, the plane pi is defined by the vertices Vi', 
V3, V4. Since the vertices involved in the edge shrinking are 
only vi, V2 if the vertices Vi, V2, V3, V4 are actually measured 
data, the sum A(pi) of the squares of the distances is 




...(8) 
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calculated by Equation ( 5 ) using an error matrix ( Kvi+Kvz ) 
relating to these vertices Vi, V2. If, in this state, the vertex 
vi' and another vertex v^, (unillustrated actually measured data) 
are converged to a vertex v/' (not shovm) to approximate the 
approximated plane pi further to an approximated plane pi' (not 
shown), the vertices involved in the edge shrinking leading to 
the approximated plane Pi' are Vi, V2, v^,. Thus, a sum A(pi.) of 
the squares of the distances is calculated by Equation (5) using 
an error matrix (Kvi+Kv2+Kv„.) relating to these vertices Vi, V2, 

Vro. 

In the second edge shrinking, it is assumed that an error 
matrix Kvi' for the vertex Vi' is (KV1+KV2), and the sum A(pi') of 
the squares of the distances is calculated by Equation (5) using 
the error matrix Kvi' for the vertex Vi' and an error matrix Kv^ 
for the vertex Vm. 

Accordingly, a general theory of the edge shrinking 
processing is as follows. If the vertices Vi, V2, V3, V4 are 
data generated by the edge shrinking instead of being actually 
measured data in FIG. 4, the error matrix Kv/ for the vertex Vi' 
defining the approximated plane Pi is Kvi ' - KV1+KV2 when the 
error matrices for the respective vertices Vi, Vz, V3, V4 are 
assumed to be Kvi, KV2, KV3, KV4. The sxim A(pi) of the squares of 
the distances is calculated by Equation ( 5 ) using an error 
matrix Kpi(= Kvi'+Kva+Kv^ = KV1+KV2+KV3+KV4 ) obtained by adding the 
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error matrices Kvi' , KV3, Kv^ for the three vertex Vi', V3, V4 
defining the approximated plane pi. 

An equation used to calculate the sum A(pi) of the squares 
of the distances is: 

KP.) = P"K,,P, •••(9) 
Here. Pi is a column vector of the approximated plane pi and mi 
is a total number of the error matrices Kv included in the error 
matrix Kpi. 

For the approximated planes p2. Pa, ...Ps. sums ACpz). 
A(P3). ... A{p8) of the squares of the distances are calculated 
by the same method as for the aforementioned approximated plane 
pi. and are substituted into Equation (8) to calculate the error 
estimation value e(Edi2) in the shrinking of the edge Edij. 

The error estimation values e(Edir). E(Ed2s). e(Edtt+i) in the 
case that the other edges Edi,. Ed^s. Ed«*i are shrunk by the same 
method are calculated. 

Subsequently, the error estimation values are sorted in 
the estimation value sorter 322 (Step #7) and it is 
discriminated whether the minimum estimation value Minjej} is 
equal to or below the permissible value inputted in Step #3 
(Step #9). If the minimum estimation value Min{eJ exceeds the 
permissible value (NO in Step #9). the data reduction processing 
by the edge shrinking is completed since a change of the shape 
resulting from application of the edge shrinking to all the 
edges is impermissibly large. On the other hand, if the minimum 
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estimation value Min{Ej is equal to or below the permissible 
value (YES in Step #9). the edge corresponding to the minimum 
estimation value Minjej} is shrunk in the edge/surface shrinking 
section 323 (Step #11). For example, if the minimum estimation 
value Min{ej} is e(Edi2) in FIG. 4. the edge Ed^ is shrunk as 
shown in FIG. 4. 

Subsequently, this routine returns to Step #7 after error 
estimation values e, (q = 1. 2. ...) are calculated for the edges 
of the polygon model after the edge shrinking that have been 
influenced by the edge shrinking (Step #13). A minimum 
estimation value Min{eq} of these estimation values is 
calculated and the edge corresponding thereto is shrunk (Steps 
#7 to #11). Thereafter, the edge shrinking is repeated edge by 
edge (loop of Steps #7 to #13), and the data reduction 
processing is completed when the minimum estimation value Min{eJ 
exceeds the permissible value (NO in Step #9). 

By performing the data reduction processing as above, the 
polygon model after the data reduction is approximated from the 
original polygon model such that the error estimation values 
from the original measurement data are equal to or below the 
permissible value. 

As described above, in the data reduction method according 
to this embodiment, the distances of the planes of the polygon 
model approximated by the edge shrinking or surface shrinking 
from the original vertices (measured three-dimensional shape 
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data) are treated as the error estimation values. Thus, the 
data reduction processing can be performed by the estimation 
standard intuitively easily understandable. 

Further, since the data reduction is performed such that 
the evaluation values lie within the specified permissible error 
range, it can be performed at a high speed while relatively 
precisely controlling the errors of the three-dimensional shape 
data forming the polygon model after the approximation. 

Particularly, reduction in the reliability of the error 
evaluation values can be suppressed even in the case that random 
noise data are mixed into the measured three-dimensional shape 
data or data reduction is performed by converging a multitude of 
vertices to one point. 

FIGS. 8 and 9 show a comparison of the error estimating 
method according to this embodiment and the conventional error 
estimating method in the case that random noise data are mixed 
into the measured three-dimensional shape data. 

FIG. 8 shows the edge shrinking of a portion (shown in 
section) where discontinuous surfaces are created on a polygon 
model by random noise data. If vertices Vi. v^ of the portion 
where large discontinuities are created on the outer surface of 
the polygon model are converged to a vertex Vi' as shown in FIG. 
8. error estimation is made based on. for example, distances 
di'. dz' of the vertex Vi' to a plane Pi4 including vertices Vi. 
V4 and a plane Pzs including vertices Vz, Vs as shown in FIG. 9 



28 




P759US 



according to the data reducing method disclosed in the above 
paper "Surface Simplification Using Quadric Error Metrics". 
Thus, the calculated error estimation values are smaller than 
sensory approximation errors. Contrary to this, according to 
the data reducing method of this embodiment, error estimation is 
based on, for example, distances di, 6.2 of the vertices Vi, Vj to 
a plane P45 including vertices Vi ' , V4 and Vg. Thus, the error 
estimation values relatively agree with the sensual 
approximation errors and are not calculated smaller than the 
latter. 

FIGS. 10, 11 and 12 show a comparison of the error 
estimating method according to this embodiment and the 
conventional error estimating method in the case that a 
multitude of vertices are converged to one point. 

FIG. 10 shows the shrinking of a zigzag surface (shown in 
section) defined by eight vertices Vi, V2, ... Vg of a polygon 
model. If the vertices Vi, V2, Ve are converged to a vertex 

Vi' located substantially in the middle between the vertices V7 
and Va as shown in FIG. 10, error estimation is made based on 
distances di', dz', ...de' of the vertex Vi' to the respective 
zigzag planes P71, P12* P23. P34* P45. Pse. Pes of the original 
polygon model as shown in FIG. 12 according to the data reducing 
method disclosed in the above paper "Surface Simplification 
Using Quadric Error Metrics". Thus, the calculated error 
estimation values are larger than sensory approximation errors. 
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Contrary to this, according to the data reducing method of this 
embodiment, error estimation is based on, for example, distances 
di, dz, ...dfi of the respective vertex Vi, V2, . . .Vg to the plane 
P78 including the vertices Vi', V7, Ve as shown in FIG. 11. 
Thus, the error estimation values relatively agree with the 
sensual approximation errors and are not calculated larger than 
the latter. 

FIG. 13 is a flowchart showing a modification of the 
procedure of the three-dimensional shape data reduction 
processing. In the procedure of the data reduction processing 
shown in FIG. 7, the estimation values e are recalculated for 
the edges influenced by the edge shrinking every time the edge 
shrinking is performed, and a next edge shrinking is performed 
based on the estimation values e. However, in the procedure of 
the data reduction processing shown in FIG. 13, the edge 
shrinking is applied to edges other than those influenced by the 
edge shrinking without calculating the estimation values e for 
the edges influenced by the edge shrinking. 

Specifically, the flowchart of FIG. 13 differs from that of 
FIG. 7 in that the operation in Step #13 is changed to Step #13a 
in which shrinking of the edges influenced by the edge shrinking 
is prohibited, and Steps #15 and #17 are added between Step #9 
and "END". 

Accordingly, common part of the flowcharts is briefly 
described and only different part is described in detail below. 
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The three-dimensional shape data of the measurement object 
G obtained by the three-dimensional shape measuring apparatus 1 
are inputted to the shape * data ' processor 3, in which a 
triangular polygon model is generated in the polygon generator 
31 (Step #1). Further, an operator inputs a permissible value 
used to judge whether the edge shrinking processing can be 
performed via the permissible value input device 33 (Step #3). 

Subsequently, estimation values £j (j = 1, 2, . . .N, N is a 
total number of edges) in the case that all edges of the polygon 
model are shrunk are calculated in the estimation value 
calculator 321 (Step #5), and the error estimation values are 
sorted in the estimation value sorter 322 (Step #7). 

Subsequently, it is discriminated whether the minimum 
estimation value Min{£j} is equal to or below the permissible 
value inputted in Step #3 (Step #9). If the minimum estimation 
value Min{Ej} is equal to or below the permissible value (YES in 
Step #9), the edge corresponding to the minimum estimation value 
MinjeJ is shrunk in the edge/surface shrinking section 323 (Step 
#11). 

Subsequently, application of the edge shrinking processing 
to the edges influenced by the previous edge shrinking 
processing is prohibited (Step #13a) and then Step #7 follows. 
This operation is made to apply a next edge shrinking processing 
to the edges not influenced by the edge shrinking processing 
this time. 
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Back in Step #7. estimation values e, (q = 1. 2. . . .Q where 
Q<N) corresponding to the edges except those for which the edge 
shrinking is prohibited are sorted, and a minimiim estimation 
value Min{£q} is calculated. It is then discriminated whether the 
minimum estimation value Min{cJ is equal to or below the 
permissible value (Step #9). If the minimum estimation value 
Min{e,} is equal to or below the permissible value (YES in Step 
#9). the edge corresponding to the minimum estimation value 
Min{e,} is shrunk in the edge/surface shrinking section 323 (Step 
#11). Further, application of the edge shrinking processing to 
the edges influenced by the previous edge shrinking processing 
is prohibited (Step #13a) and then Step #7 follows. 

Thereafter, by repeating the above operations, a next edge 
shrinking processing is performed while excluding the edges 
already influenced by the previous edge shrinking processing 
every time the edge shrinking processing is performed (loop of 
Steps #7. #9. #11, #13a). 

On the other hand, if the minimum estimation value Min{eq} 
becomes larger than the permissible value (NO in Step #9) while 
the edge shrinking processing is repeated, it is discriminated 
whether the edge shrinking has been prohibited for these edges 
(Step #15). If the discrimination result is affirmative (YES in 
Step #15). Step #7 follows after error estimation values (r = 

2. ) are recalculated for these edges (Step #17). A 

minimum estimation value Min{eJ of the error estimation values 
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ar0 calculated, and the edge corresponding to the minimum 
estimation value Min{6r} is shrunk (Steps #7, #9. #11). 

If the minimum estimation value Min{eq} becomes larger than 
the permissible value and there is no more edge for which the 
edge shrinking is prohibited (NO in Steps #9. #15). the data 
reduction processing by the edge shrinking is completed. 

FIG. 14 shows an exemplary data reduction process by the 
operation procedure shown in FIG. 13. In FIG. 14. a polygon 
model Mo is an original polygon model. This polygon model Mo is 
formed by arraying 28 vertices Vy. v^. . . .v^s in 7 rows such that 
the numbers of the vertices in the respective rows are (3. 4. 5. 
4. 5. 4, 3) and connecting adjacent vertices. 

Edges el. e2. e3. e4. e5 are selected as such to be 
successively shrunk by the operation procedure shown in FIG. 13 
to simplify the description since they are not mutually 
influenced in the case that the respective edges ei (i = 1. 
2 . . . . 5 ) are shrunk . 

Accordingly, the data reduction process of FIG. 14 shows a 
case where the edges el. e2. . . .e5 of the polygon model Mo are 
successively shrunk by the operation procedure of FIG. 13. 
Polygon models Mi. M^. M3. M4. M5 are approximated polygon models 
obtained when the polygon model Mo is successively approximated 
by successively shrinking the edges el. e2. ...e5. 

According to the procedure of the data reduction processing 
shown in FIG. 13. when the polygon model Mo is first 
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approximated to the polygon model Mi by shrinking the edge el by 
the first edge shrinking processing, edges formed by a converged 
vertex Va and its neighboring vertices Vg, Vio. Vn, V13, Vie, Via, 
Vi9, V20 {vertices located within a dotted- line frame in the 
polygon model Mi) are influenced. Thus, in the next edge 
shrinking processing, the estimation values Cq are sorted again 
for the edges (those located outside the dotted-line frame of 
the polygon model Mi) excluding the above edges, and the polygon 
model Ml is approximated to the polygon model M2 by shrinking the 
edge e2 corresponding to the minimxim estimation value Min{Eq}. 

Further, since the edges defined by a converged vertex Vb 
and the vertices Vi, Vz, V4, Vg, Vb, Vio, V13, Va (vertices located 
within a dotted-line frcune in the polygon model M2) are 
influenced in the next edge shrinking processing, the estimation 
values Es are sorted again for the edges (those located outside 
the dotted-line frame of the polygon model M2) excluding the 
above edges , and the polygon model M2 is approximated to the 
polygon model M3 by shrinking the edge e3 corresponding to the 
minimum estimation value Min{es} . 

The edge shrinking is applied to the edges e4 , e5 by the 
same method, and the original polygon model Mo is approximated 
to the polygon model M5 and the number of data is reduced by 5 
by the five edge shrinking processings. 

If a process for approximating the original polygon model 
Mo to a 6*'*' polygon model Mg (not shown) is considered, no 
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estimation value is recalculated for the edges having been 
influenced by the edge shrinking according to the operation 
procedure shown in FIG. 13. Thus,* in the example of FIG- 14, 
when the estimation value is first calculated for 66 edges 
forming the original polygon model Mo, the evaluation values are 
recalculated for 50 edges forming the polygon model Mg without 
calculating the estimation values again for the approximated 
polygon models Mi to M4. Therefore, the estimation values are 
recalculated 50 times until the 6^^ polygon model Ms is obtained. 

On the other hand, since the estimation values are 
recalculating for the edges influenced by the edge shrinking 
according to the operation procedure shown in FIG. 7, the 
estimation values are recalculated for 16 edges influenced by 
the edge shrinking for each of the approximated polygon models 
Ml to M5. Thus, according to the operation procedure shown in 
FIG. 7, the estimation values are recalculated 80 (16x5) times 
until the 6*^ polygon model Me is obtained. Therefore, the 
operation procedure shown in FIG. 13 has an advantage of high- 
speed processing since the estimation values are recalculated a 
reduced number of times . 

FIG. 15 shows a modification of the flowchart of FIG. 13. 
In the flowchart shown in FIG. 15, the edges influenced by the 
edge shrinking are made not to be shrunk in the next edge 
shrinking by replacing the estimation values thereof by a 
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predetermined estimation value larger than the permissible 

value . 

Specifically, the flowchatt of FIG. 15 differs from that of 
FIG. 13 in that Step #13a is changed to Step #13b in which the 
estimation values of the edges influenced by the edge shrinking 
are set at £„^ax. and the discrimination criteria in Step #15 and 
the objects in Step #17 whose estimation values are recalculated 
are changed to "edges whose estimation values are Cmax" . and Steps 
#19 and #21 are added after Step #17a to shrink the edges whose 
estimation values are recalculated. Accordingly, only part of 
the flowchart of FIG. 15 different from that of FIG. 13 is 
described in detail below. 

When the edge corresponding to the minimum estimation value 
Min{Ej} is shrunk in the edge/surface shrinking section 323 (Step 
#11), the estimation values of the edges influenced by this edge 
shrinking processing are set at the predetermined estimation 
value e^ax larger than the permissible value (Step #13b) and then 
Step #7 follows. This operation is performed so as not to 
shrink these edges in the next shrinking processing , i.e., to 
apply the edge shrinking only to the edges not influenced by the 
edge shrinking processing this time. 

Back in Step #7 , evaluation values e± (q = 1 * 2 , — N-1 ) 
are sorted, and a minimiam estimation value Min{£i} is calculated. 
It is then discriminated whether the minimum estimation value 
Min{ei} is equal to or below the permissible value (Step #9). If 
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the minimum estimation value Min{ei} is equal to or below the 
permissible value (YES in Step #9), the edge corresponding to 
the minimum estimation value Min{ei} is shrunk in the 
edge/surface shrinking section 323 (Step #11). 

Subsequently, the estimation values of the edges influenced 
by the edge shrinking processing this time are set at Emax (Step 
#13b) and then Step #7 follows again. Thereafter, by repeating 
the above operations, the next edge shrinking processing is 
performed while setting the estimation values of the edges 
influenced by the edge shrinking process at Emax every time the 
edge shrinking processing is performed (loop of Steps #7, #9, 
#11, #13). 

On the other hand, if the minimum estimation value Min{ei} 
becomes smaller than the permissible value (NO in Step #9) while 
the edge shrinking processing is repeated, it is discriminated 
whether there is any edge having the estimation value of Emax 
(Step #15a)- If there are such edges (YES in Step #15a) , the 
estimation values Es (s = 1, 2, ...) are calculated for these 
edges (Step #17a). 

A minimiim estimation value Min{£s} is calculated by 
discriminating the error estimation values Es (Step #19). If the 
minimum estimation value Min{Es} is equal to or below than the 
permissible value (YES in Step #21), Step #11 follows in which 
the edge corresponding to the minimum estimation value Min{£s} is 
shrunk . 
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On the other hand, if there is no more edge having the 
estimation value of e,^ (NO in Step #15a) or the minimum 
estimation value Min{es} of the estimation values Es calculated 
for the edges having the estimation value of becomes larger 
than the permissible value (NO in Step #21) while the edge 
shrinking processing is repeated, the data reduction processing 
by the edge shrinking is completed. 

Since the shrinking of the edges influenced by the edge 
shrinking is prohibited in the operation procedure of FIG. 15 as 
in the operation procedure of FIG. 13 by setting the estimation 
values of such edges at the estimation value E,„ax larger than the 
permissible value, the estimation values are recalculated a 
reduced number of times and high-speed processing can be 
performed. 

In the procedures of the data reduction processing shown in 
FIGS. 13 and 15. the estimation values are not calculated for 
the edges influenced by the edge shrinking processing, i.e.. 
these edges are not shrunk until the edge shrinking processing 
is prohibited for all edges. However, the prohibition of the 
edge shrinking processing may be released if a specified 
condition is met such as when the edge shrinking processing is 
repeated a specified number of times or when a difference 
between the minimum estimation value and the permissible value 
lies within a specified range, and the edge shrinking may be 
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made possible again by recalculating the estimation values for 
all edges. 

Such an arrangement enables the data reduction processing 
to be applied to the edges for which the edge shrinking 
processing is prohibited, but the data reduction processing can 
be performed within a permissible error range. Therefore, data 
can be maximally efficiently reduced within a permissible range. 

Since the procedures of the data reduction processings 
shown in FIGS. 13 and 15 are designed to increase the processing 
speed by reducing the number of times the estimation values are 
recalculated, the estimation values are not limited to the 
errors £ expressed by Equation (6). The data reduction 
according to the present invention is widely applicable even in 
the case that estimation values are changed by the edge 
shrinking or surface shrinking and need to be recalculated after 
shrinking, for example, those disclosed in the aforementioned 
paper "Surface Simplification Using Quadric Error Metrics" are 
adopted . 

Although reduction of the three-dimensional measurement 
data obtained by the three-dimensional shape measuring apparatus 
1 is described in the foregoing embodiment, the present 
invention is not limited thereto and is widely applicable to 
three-dimensional data created by computer graphics or the like. 

As described above, the three-dimensional shape data are 
reduced by converging two or more vertices to one vertex to 
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shrink the edges or surfaces. Specified estimation values are 
calculated for all the surfaces deformed by shrinking the edges 
or surfaces of the polygon model based on the distances between 
the respective surfaces and all the original vertices involved 
in the deformation of these surfaces. The three-dimensional 
shape data are reduced by edge shrinking or surface shrinking 
when the estimation values are equal to or below a predetermined 
permissible value. Accordingly, the error can be controlled by 
the estimation standard intuitively easily understandable, i.e., 
by approximation errors from the original three-dimensional 
shape data. 

Also, even in the case that, for example, the orientation 
of the surface locally experiences a sudden change into a zigzag 
surface in the polygon model, the data reduction processing can 
be suitably performed while the errors are relatively precisely 
controlled. 

Particularly, since the estimation values are assumed as 
errors e defined by Equation (6), the errors e can be easily 
calculated only by storing an error matrix calculated for the 
original three-dimensional shape data. Thus, a high-speed 
processing is possible and a memory necessary for the 
calculation of the data reduction processing can be efficiently 
utilized. 

Furthermore, the three-dimensional shape data is reduced by 
repeating the edge or surface based on estimation values. 
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However, the application of the shrinking to the edges 
influenced by the previous edge shrinking is prohibited. 
Accordingly, recalculation of estimation values can be executed 
a reduced number of times, which thus enables a higher speed 
processing. 

The converged vertex portions, which are influenced by the 
edge or surface shrinking, are given with a predetermined value 
to place these portions out of the shrinking, which can thus 
assure easier shrinking prohibition to a particular portion. 

Moreover, in the case that a specified condition is 
satisfied, the prohibition is released to allow the shrinking 
processing to the prohibited portion again. Accordingly, data 
reduction can be accomplished as efficiently as possible within 
a permissible range. 

As this invention may be embodied in several forms without 
departing from the spirit of essential characteristics thereof, 
the present embodiment is therefore illustrative and not 
restrictive, since the scope of the invention is defined by the 
appended claims rather than by the description preceding them, 
and all changes that fall within metes and bounds of the claims, 
or equivalence of such metes and bounds are therefore intended 
to be embraced by the claims. 
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